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Abstract 

Inspired by the discrete evolution implied by the recent work on loop quantum cosmology, we 
obtain a discrete time description of usual quantum mechanics viewing it as a constrained system. 
This description, obtained without any approximation or explicit discretization, mimics features of 
the discrete time evolution of loop quantum cosmology. We discuss the continuum limit, physical 
inner product and matrix elements of physical observables to bring out various issues regarding 
viability of a discrete evolution. We also point out how a continuous time could emerge without 
appealing to any continuum limit. 
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I. INTRODUCTION 



The recent work of Bojowald on Loop Quantum Cosmology (LQC) obtains the Wheeler- 
DeWitt equation (Hamiltonian constraint) as a difference equation [|l|. The eigenvalues 
of the volume operator are discrete in quantum geometry and are taken as playing the 
role of a 'discrete time' in the context of isotropic LQC [Q]. The order of the difference 
equation is typically high (16 for isotropic, Bianchi-I) and the number of independent (and 
non-degenerate) solutions is reduced by one due to the coefficient of the highest (lowest) 
order term vanishing for discrete time equal to zero [Q. This feature is crucial for the 
'singularity avoidance' mechanism [Q. Furthermore, a continuum limit is defined wherein 
the Immirzi parameter plays a crucial role. This limit is used to distinguish the so called 
'pre-classical' solutions and it is shown that the singularity avoidance mechanism also leads 
to a unique (up to normalization) pre-classical solution of the Wheeler-DeWitt equation 
P, ^. These are very interesting indications that the procedures of loop quantum gravity 
as applied to cosmological mini-superspaces, do lead to physically reasonable solutions. 
Although certain choices of definitions of constraint operators with certain factor orderings 
can be made with reasonable justification, the procedure is not devoid of ambiguities 0. 
While there are many issues to be resolved, we focus on one feature, namely 'discrete 
evolution', which seems to be robust. 

For instance, mathematically, solutions of the Wheeler-DeWitt equation can be pre- 
sented as a sequence of states (isotropic LQC). The sequence label is highly suggestive 
of a discrete 'time'. This dynamical interpretation is something additionally attempted 
and its viability needs to be established. In a quantum theory, an evolution interpretation 
must be established at least at the level of expectation values, the states being not directly 
observable. Furthermore it is not enough to generate a family, continuous or discrete, of 
expectation values. It should be possible to detect the changes paying due attention to 
the uncertainties. By contrast, a continuum limit does not appear to be essential for a 
dynamical interpretation even though emergence of a continuous time description in a 
semi-classical limit is of course desirable. Presently, in the context of LQC, these issues are 
discussed somewhat schematically. One of the motivations for the present work is to have 
simple, well known examples which nonetheless mimic steps taken in LQC. 
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It is well known that the usual quantum mechanical systems (or classical for that matter) 
can also be viewed as a constrained system leading to a frozen-time description. The Dirac 
procedure nonetheless allows one to interpret the physical state condition as an evolution 
equation - the usual Schrodinger equation. By using the 'Fock' representation instead of 
the Schrodinger representation for the extra time degree of freedom, one can get a discrete 
'evolution' equation which mimics all the features seen in the LQC work. The order of the 
difference equation is two but there is a reduction of number of non-trivial independent 
solutions; there is a parameter analogous to the Immirzi parameter which can play a similar 
role in exploring pre-classicality and a continuum limit. 

To explore evolution at the level of expectation values of course needs definitions of 
physical inner products, observables and their matrix elements. The advantage in the 
quantum mechanical case however is that we know physical inner products and physical 
matrix elements so that we can push further the interpretation of the difference equation 
as "really" an evolution equation. Furthermore since the Schrodinger and the 'Fock' 
representations are equivalent, we can relate the continuous and the discrete evolutions by 
a transform. It turns out that the dynamical interpretation is not as straightforward as 
indicated by the LQC works. 

'Discrete Time' has appeared in the literature several times in various forms and 
with various motivations. The present work is very different from these earlier works. In 
particular, we are not seeking a discrete time formulation, ab initio, for any particular 
reason. We observe that in a frozen time formulation of dynamics, a 'time' appears as a 
basis label which has an arbitrariness about it. The dynamics is then obtained as a family 
of states labeled by the 'time'. One has natural choices of continuous and discrete labels. 
The choice of the continuous label leads to the usual quantum dynamics and we explore 
the discrete choice in detail, in particular with regards to observability of evolution. 

In section II we detail the case of usual quantum mechanics cast in a frozen time form, 
both with continuous and discrete time and exhibit its analogy with LQC. We discuss 
the continuum limit and show the relation between the continuous and the discrete time 
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descriptions. In section III we discuss natural candidates for physical quantities needed 
to push the evolution interpretation further and the difficulties encountered. In the last 
section we discuss the issue of interpretation in some generality and point out a possible role 
of the parameter (3 appearing in the discrete description. We conclude by making a series 
of remarks. 

II. NON-RELATIVISTIC QUANTUM MECHANICS 
A. Frozen Time Description: Continuum Case 

Let To denote a classical phase space and let Tio denote the Hilbert space of the corre- 
sponding quantum system. Let r, vr,- denote two extra phase space variables corresponding 
to the usual continuous time and its 'conjugate momentum'. Let '■= x Lq denote 
an extended phase space and let Hkin '■= L2{R,dT) ® Ho the corresponding kinematical 
Hilbert space. At the classical level we impose the single constraint := vr^ + H{uj) ~ 0. 
Here, u denote the usual phase space coordinates and H{uj) denotes a time independent 
Hamiltonian on Tq. Quantum mechanically, the operator version of the constraint is 
imposed on to select the 'physical states'. Explicitly, the physical states are those on which 
the operator tTt- + H vanishes. The constraint operator is of course is not identically zero. 

2 

For definiteness one may think of Tq = R'^'^ ,Ho = L'^{q\d"'q) and H^u) = ^ + V{q^) 
though this is not necessary. There is no external time any more. 

At the classical level, the Dirac observables, y4(r, tTt-, cj) are defined by {A,(j)}pB ~ 0. 
Some simple examples of Dirac observables are: functions of only vr^ and functions indepen- 
dent of r, tTt- which Poisson commute with the Hamiltonian (in particular the Hamiltonian 
itself). This is a rather limited class of observables. One could however choose a r dependent 
family of functions of lu satisfying the differential equation |^ + {A{t,uj), H{uj)}pb = 
with the initial condition A{0,uj) = Aq{uj). Such solutions of the differential equation are 
trivially Dirac observables. In particular, the usual solutions of Hamilton's equation are 
also Dirac observables but with r —r. These are the 'evolving' observables in a frozen 
time formulation |0. 



4 



In the quantum description, one quantizes the r, vr,- in the usual manner and the constraint 
equation becomes just the Schrodinger equation. Exphcitly, one can write a general vector 
in Tikin in the form: 

I*) = J dr\T) ® |$(r)) (1) 
The kinematical inner product is then given by, 

(^'|^)..n = I dT{<i>\r)mr))o (2) 
The suffix refers to the inner product in TYq- 

Thus physical states are those whose |$(t)) satisfy the usual time dependent 
Schrodinger equation. Introducing f/(r) := e"'^'^^, we denote the solutions of the 
Schrodinger equation as |$(t)) = f/(r)|$(0)). The corresponding |\l/)'s are not normaliz- 
able with respect to the kinematical inner product since the integrand is independent of r 
rendering r integration divergent. 

Dirac observables are usually defined as those observables which commute with the con- 
straints. It may be sufficient to require that the physical observables commute only weakly 
with the constraints. In this case the factor ordering must ensure that the constraint opera- 
tors act first i.e. are to the right. Weak commutation then amounts to physical observables 
acting invariantly on the physical states. As an example one can define 'evolving observables' 
as follows 0. Corresponding to a usual operator, O on Hq, define a family of operators on 
the space of physical states by, 

[d{T)\mr'):=U{T')U{-T)dmT) (3) 

These 'evolving observables' are physical in the sense they act invariantly on the space 
of physical states. In general, they do not commute with the constraint operator in the full 
Ti-kin- Since physical states are not kinematically normalizable, one has to define a physical 
inner product on the space of solutions of the constraint. A natural definition suggests itself 
as: 

:= (<I>'(ro)|<I>(ro))o tq is some fixed time (4) 

The kinematical inner product is just the integral of the physical inner product over tq. 
It is obvious that the physical inner product is independent of the particular tq chosen and 
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of course coincides with the usual inner product in Ho- 



lt is easy to see that the physical matrix elements of these evolving observables between 
physical states are given by, 

{^'\d{T)\^Uy = {^'mu{-T)du{T)\m)o (5) 

Hence, the physical matrix elements are obtained as the usual matrix elements of oper- 
ators on Hq. We see thus that the usual description of quantum dynamics in terms of a 
continuous time can be recast as a frozen time presentation. This is of course well known. 
We will now introduce a discrete time description which mimics all the features seen in the 
loop quantum cosmology. 

B. Discrete case 

Instead of choosing the usual Schrodinger representation for r, tt^ = l^t us intro- 

duce a number representation. Define a :— aT+i^TTr and its hermitian conjugate, a, (5 are 
real and satisfy 2aj3h = 1. Hence we have a one parameter family of creation-annihilation 
operators labeled by (3, say. This parameter is expected to play a role analogous to that 
played by the Immirzi parameter in LQC. We will choose the eigenvalues of the number 
operator, N :— a^a as our discrete time label. Notice that these eigenvalues are independent 
of p. For future use in continuum limit we note that. 

This will justify the eigenvalues of the number operator being identified with r at least 
for large eigenvalues and small P and of course n is monotonic in r^. 

In terms of creation-annihilation operators, the constraint equation becomes. 
Writing |^) = Yl,^=o 1^) ® |$n) one obtains the kinematical inner product as, 

oo 
n=0 
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The physical states are then those |^) whose |$n)'s satisfy, 

This is our discrete time evolution (discrete Schrodinger equation). We notice that 
the difference equation is an operator difference equation of order two implying that two 
vectors (in Ho ) have to be specified to determine a solution i.e. has 'two independent' 
solutions. However exactly as in the case of loop quantum cosmology, we have a consistency 
condition (since spectrum of N is bounded below) which fixes |$i) in terms of |$o) 
implying a unique solution for every given |$o)- By contrast, in the continuum descrip- 
tion, there is no such condition but the equation is of course a first order differential equation. 



For subsequent analysis, it is convenient to convert the vector equation in to infinitely 
many scalar equations by expanding the = ^n\^a) in the eigenbasis of the H. For 
simplicity we have assumed that spectrum of H in TIq is discrete. The discrete equation 
then becomes. 



- -tSi^-' W^^" V n > -1, V „. (10) 

Let us now turn to the "continuum limit". This can be understood in various ways. A 
simple way is to ask whether one can find continuous variable t and a function C"(t) which 
will interpolate a solution of the discrete equation for large n. To explore this let us look for 
a function C{t) and a function t{n) such that in a suitable large n limit one has (suppressing 
the label a), 

dC 

Cn+k-=C{t{n + k)) := C{t{n) + kSt) ^ C{t) + kSt—. (11) 

where, kSt — t{n + k) — t{n) has been used. Treating 5t small is equivalent to requiring 
t{n) to be slowly varying with n. Substitution in the difference equation and keeping terms 
to leading order in n (treating 6t also as small), leads to C^^^ = — lEj^j^. The left 
hand side is a function of t by assumption so we must have j^j^ to be a function of t. Since 
the limit is to be considered for all C", we have excluded £. Simplicity and dimensional 
considerations then suggest that we choose t such that :— h~^. The variable t so 
specified will be denoted by r. It follows that T(n) = 2hl3y/n. Observe that 5t ~ rT^^'^ and 
thus vanishes for large n without having to take h(5 to be vanishingly small. This is different 



from the LQC To get a finite r for arbitrarily large n however, we must consider a joint 
limit n —>■ oo,/3 —>■ keeping r fixed. We could have taken hjS to zero but presently we are 
interested in a continuum limit instead of a semi-classical limit. It follows that in the above 
joint limit C(r) satisfies the usual Schrodinger equation and so does the full |$)(r). Thus 
we see that continuous functions that interpolate solutions of the difference equation asymp- 
totically are solutions of the usual Schrodinger equation with a suitable identification of T{n). 

Can we obtain C(r) as a limiting function from the joint limit of C„? The answer is yes. 
Consider the difference equation for large n but with fixed f3, 



The left hand side of this equation equated to zero is an equation with constant 
coefficients (Poincare type) . The asymptotic behaviour of its solutions is given in terms 
of its characteristic roots obtained by substituting ~ A". The characteristic roots are 
just A = ±1, independent of the label a. Evidently, the root A = — 1 can not correspond 
to a solution which has limiting value in the joint limit. Furthermore, even for A = 1, one 
can not see the limit to be a solution of the continuum Schrodinger equation. One needs a 
more refined ansatz: C„ ~ A"/i^. Substitution determines A = ±1 from the n independent 
term and also gives infx = —2i\(3Sa from the sub-leading n"^/^ term. For A = 1 we see that 
goes over to the solution of the continuum Schrodinger equation. By contrast , A = — 1 
does not have a limit. A generic asymptotic solution will be a linear combination of these 
two asymptotic solutions and it will not have a limiting value in the joint limit. There is 
then a unique solution that does have limiting value. This is very similar to the arguments . 

So far the steps are completely analogous to those taken in loop quantum cosmology. 
However here we run in to a problem. No exact solution of the difference equation can 
possibly have a non-zero and finite limiting value in the joint limit. This follows because 
the ratio of the for n odd and n even is necessarily purely imaginary as is evident 
from the difference equation. Thus although the assumptions we made about a conceivable 
continuum limit do admit a corresponding ansatz for the asymptotic solution of the 
difference equation, no exact solution can in fact support such an ansatz. The notion of 
pre-classical limit as articulated in LQC 0, is not realized by any exact solution even 




(12) 
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though a continuum hmit in the sense of difference equation going over to a differential 
equation is vahd. One has only asymptotic solutions i.e. solutions of asymptotic equation as 
distinct from asymptotic form of solutions of the exact equation, which have a pre-classical 
limit but no exact solution has this property. It seems that the continuum solutions can at 
best be thought of as approximating the exact solution that too only asymptotically i.e. for 
large n. 



Is not having any pre-classical solutions a disaster for getting a continuum picture? Not 
necessarily. For emergence of a continuum description from an underlying discrete one what 

is needed is a mapping to continuous description and not necessarily a continuum limit. The 
next sub-section shows how this can happen. 



C. Relating the continuum and the discrete descriptions 

The evolution equations have been derived by writing: 

POO 

1^) e Hkin — / dr\r) |$(t)) continuous case 

J oo 

oo 

= ^ \n) (g) \^n) discrete case (13) 



n=0 



Imposing the constraint in the Schrodinger representation and the Fock representation 
respectively leads to the usual Schrodinger equation for |$(t)) and the difference equation 
for |$n). The two basis vectors \r) and \n) are related as: 

oo 

|r) = El^)(^l^) = E/n(r)|n), 

n=0 n 

\n) = r dr\r){r\n)^ [ drfnir)\r) (14) 

J — oo J 

The transformation functions, /n(T) are easily determined and are given by. 



a/2 hf3'> 



Nn = ((27r)V4y^v^) ' (15) 

and i?n(0 the Hermite polynomials. The vectors |$(t)) and are easily seen to 
be related as: 
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Mr)) = 5^/„(r)|<|.„) 

n 

l$n) = I rfr/:(r)|$(r)). (16) 

Using the properties of the Hermite polynomials, it is easy to see that |$„,) satisfy the 
difference equation Zj^|$(r)) satisfies the Schrodinger equation. The respective initial states 



|y^"rfee-^-^^«^^||$(0)). (17) 



are related as: 

By expanding |$(0)) in eigenstates of the Hamiltonian one can show that, 

($o|$o)o < 2v^/i/3($(0)|$(0))„. (18) 
This shows that in the limit /9 — > the transform breaks down. 

This route, available in the present case, shows that it is possible to define states depend- 
ing on a continuous variable in a differentiable manner without appealing to any 'pre-classical' 
or otherwise limiting procedure. 

III. PHYSICAL QUANTITIES 

As a first attempt, we will just mimic the steps followed in the continuous time case. 
For this it is convenient to write the second order difference equation as a first order matrix 
difference equation. This is easily achieved [H. The evolution equation can be written as: 



Zn+i = A{n)zn n > 
z^_^ = M{n)zn n>l. (19) 



It follows from the definitions that A{n — l)M{n) = M{n + l)A{n) = I. Not all states 
evolving by the above equations are physical though because the physical states also have 
to satisfy the consistency condition namely |$i) = —2i(3H\^Q). 



1 l^n+l) 1 








1 |$n-l) ^ 
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Let is denote by D„ the space spanned by the z^s. It can be viewed as Ho ^Ti-o- To 
define an evolving observable 0(m) corresponding to an operator O acting invariantly on 
all DnS, consider the Zm member of an physical state {|<l>„)}. Operate on it by Oq and 
evolve back to Dq. This will not in general satisfy the consistency condition. Let V be 
a projection operator which will project any element of Dq on to one corresponding to a 
physical one. That is, define a new zq whose components satisfy the consistency condition. 
Evolve this to any n > 0. This clearly defines a physical state. 0{m) is defined to produce 
this physical state from the starting physical state. By construction, these operators, defined 
to act on physical states, produces a new physical state. Introduce the evolution operators 
E{0, n) : Dq —>■ D^ and E{n, 0) : D^ Dq given explicitly by, 

E{0,n) := A(n- l)A(n-2)---A(0) 

E{n,Q) := B(1)B(2) ■ ■ -1(71) 

Then 0{m) is expressible as, 

[d{m)z]n := E{0, n)oVo E{m, 0)dzm 
The projection operator can be constructed easily and is uniquely given by. 

In analogy with the continuum case, one may naturally define a physical inner product as 
the n = term of the series in equation 8. We can make other possible choices, for example, 
z^Mzq, where is a suitable 2x2 matrix of operators on Hq. M. = diag{I, O) will produce 
the n = term of the series. Physical matrix elements of the evolving observables are then 
given by, 

= ZQ^VMVE{m, 0)dE{0, m)zo (23) 

We have used ZqV = Zq for physical states to obtain a symmetrical expression. While 
this looks very similar to the continuum case (apart from the presence of VAiV), its 
implications are very different. 

These definitions, while plausible, are unsatisfactory. We have focussed on the operators 
which act invariantly on the space of physical states. However to be of use for measure- 
ments, such operators must satisfy further properties such as self-adjointness. This of 
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(20) 
(21) 



course needs to be defined relative to the physical inner product. The non-unitarity of 
evolution however implies that even if have a self-adjoint operator at n — 0, the other 
members of the family are not self-adjoint in general. The presence of the projection 
operator implies that algebraic relations among operators are not preserved by the 
evolution. For example, an operator associated with an 0^(0) is not the square of the 
corresponding operator associated with 0(0) and like wise for commutation relations. 
In the present case of Fock representation, the condition arises due to the spectrum 
of number operator being bounded below. In LQC, although the state label n takes all 
integral values, there is still the consistency condition which should cause similar difficulties. 

One could have dealt with the second order equation it self. Now there is no need for 
any explicit projection operator. The evolution is still non-unitary but in addition, one 
can not evolve a given |$„) back to a |$o) since the equation is second order. The first 
order formulation avoids this but introduces explicit projection operator. Thus our attempt 
to mimic the steps followed in the continuous time case do not lead to satisfactory definitions. 

However, wc can appeal to the relations between the continuum and the discrete descrip- 
tion discussed before. Then the physical inner product and matrix elements as defined in 
the continuous case can be expressed as, 

(*'(r)|$(r)) = ^/;(r)/„(r)(*:„|*„) 

m,n 

($'(r)|/l|$(r)) = E/4(r)/.(r)(Cl^l*n) (24) 

These are very different from the physical inner products and matrix elements we at- 
tempted previously! The right hand sides involve infinite sums and are highly non-local in 
the discrete time label. Note that the apparent r dependence on the right hand sides is 
consistent with that implied by the left hand sides of the above equations. Thus, if we 
somehow invented these inner products and matrix elements, we could construct a contin- 
uum description from the discrete one. In the case of quantum mechanical example we are 
discussing, we have the advantage of knowing a continuum description ab initio but for LQC 
also something similar can be conceivable. This however is not attempted in the present 
work. In the next section we make some general remarks regarding viability of a 'dynamical' 
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interpretation and attempt to arrive at an interpretation of (3. 



IV. 'EVOLUTION' IN QUANTUM MECHANICS 

Unlike classical mechanics, quantum mechanics permits two notions of evolution which 
are not equivalent due to the uncertainty relation. The two notions correspond to evolution 
at the level of states i.e. a continuous or a discrete family of rays (or vectors) and 
evolution at the level of observed quantities i.e. family of expectation values of observables. 
To appreciate the non-equivalence of these two let us quickly recall the derivation of 
time-energy uncertainty relation. One can do this quite generally. 



Let G be a self adjoint operator on a Hilbcrt space. Consider the one parameter group of 
unitary operators generated by G, U{^) := exp{—i^G) , ^ eM.. Define a family of normalized 
state vectors ^(OIV'lO))- For ^-^y self adjoint operator, A, corresponding to an 

observable define '■— ('0(OI^IV'(O)- Assume A to be independent of ^ for simphcity. 

Then it follows that, 

5U{0:^SC^^i5mO\[G,A]m)) , > (say) (25) 

Defining G' :— G — ('0(OI^|V'(O) like wise for A one gets, 

\8fm = \ m{mG',mm i 

< 25i |Im(G'V,^V)| 

< 25^ |(GV,^V)I 

< 25i\\G'm\ m 



< 25C AG^ AA^ (26) 

Clearly in order to detect a change in the expectation value fipiC), the change computed 
above must be at least as large as the uncertainty, AA^. This immediately gives the uncer- 
tainty relation: 

AG^ > ^ (27) 

Note that this derivation is independent of canonical commutation relations and 
thus is not tied to a phase space ~ M^^, 6^ ^ A^^ although it could be. The above 
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derivation has also assumed 5$, to be small enough so that the higher order terms can be ne- 
glected. If these are also included, then the uncertainty relation will assume a different form. 

Taking G to be the Hamiltonian and ^ = r/h we get the usual time-energy uncertainty 
relation while for G equal to the momentum (say) and ^ = q/h we get the position- 
momentum uncertainty relation. Its meaning is that we can not observationally resolve 
values between ^ and ^ -|- 5^. The fact that there is a non-zero lower bound and AG^, in 
a physical situation is always finite (though it could be very large) implies that continuum 
values for ^ are strictly mathematical idealizations. This distinguishes the two notions of 
'evolution' in quantum mechanics, mentioned above. The states can be thought of as evolv- 
ing continuously but observationally, continuous evolution is necessarily an idealization. 
The absence of a non-zero lower bound in classical mechanics permits continuum val- 
ues of ^ to be taken more hterally. Note that this applies not just to 'time' but also to 'space'. 

The notion of observationally detectable evolution can be articulated as follows. We 
can meaningfully say that a system has changed its state provided we can measure at least 
one of its properties and detect a change. Any such measurements will give expectation 
values together with uncertainties for the corresponding observable. Thus to conclude that 
a system in some given state has changed 'over a period of time' one must be able to 
find at least one observable whose expectation value in that state changes more that the 
uncertainty. Note that this must be understood at the level of an ensemble of identically 
prepared systems since a single measurement on a single system will just produce some 
eigenvalue of the observable according to the standard interpretation of quantum mechanics. 
To account for 'over a period of time', one must assume a family (discrete or continuous) 
of states in which the expectation values are to computed. 

Thus, for the observational notion of an evolution, the central quantities are expectation 
values. Given a discrete family of vectors and self adjoint operators (or a family thereof) 
one can construct a corresponding family of expectation values. Such a family could be 
usefully interpreted as an 'evolution' provided that the difference between consecutive 
members of the family of expectation values is larger than the corresponding uncertainties, 
at least for some observable and for generic states. Such a criterion of detectable evolution is 
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independent of how the family of vectors is chosen and whether the members of this family 
are connected by unitary operators. This is also independent of whether the sequence of 
vectors is obtained from solutions of some (Hamiltonian) constraint of a constrained system. 



As seen in the example of quantum mechanical system, the net result of imposing con- 
straint in the kinematical arena is to produce a family of vectors in TIq. Since such families 
are uniquely determined by |$o) (or 1^(0))), the space of 'physical states' is isomorphic to 
TYo- A natural choice of physical inner product is then simply the inner product in TYo and 
'physical observables' are naturally self adjoint operators on this Hilbert space. We can now 
construct a sequence of 'physical" expectation values: 

Expanding the states in terms of eigenstates of the Hamiltonian (assuming discrete spec- 
trum for simplicity), = Yl,pCn,p\^p)-, we get, 

{A)„ = S>.-^j;^^;:;'<^;l-^l^-> (29) 

Z^p \^n,p\ 

If |$o) is an eigenstate of the Hamiltonian, then so are |$„) V n. The expectation 
values defined above are then independent of n. Thus eigenstates of the Hamiltonian are 

'stationary' states even with respect to the discrete 'evolution'. 

The coefficients Cn,p above, satisfy the difference equations and are n*^ order polynomials 
in {—2i(5Ep). The expectation values and the uncertainties are thus rational functions of (5. 

Applying the reasoning to the simplest two level system brings out further possibilities. 
Let TYo be two dimensional and let the Hamiltonian he H = £a^. Let a generic observable 
be a Hermitian 2x2 matrix. Then it is easy to see that all expectation values (and hence 
also the uncertainties) are independent of f3\ Explicitly, let 

(cos9 \ I a b \ 

A -.^ \ a, c real. (30) 

sin9e''t' I \b*cj 

Putting |$„) := -^|$o), Vn > and Q :— —2il3H, one can see that the P„'s satisfy the 



equation, 

Pn+2 = QPn+1 + in+ l)Pn, Vn > -1, Pq := I, P-1 := (31) 
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By inspection, P„'s are n*'* order polynomials in Q. Furthermore, for even (odd) n, only 
even (odd) powers of Q occur. Since is a multiple of the identity matrix, is a 
polynomial in (— 4/9^£^) times the identity matrix while P2m+i is another polynomial in the 

same variable times (73. In computing the expectation values, these polynomials cancel out 
in the ratios and all the j3 dependence disappears leaving us with, 

{A)n = a cos^e + c sin^e + {-Ifsine cosO {be''^ + 6*6"^'^) 
(A^)^ = \bf + siri'e cos^e {{a - cf - (6e'^ + b*e-"t'f} 

+ {-lysine cose cos2e (c - a)(6e^'^ + 6*6"^'^) (32) 

Prom these it follows that change in the consecutive expectation values depends only on b 
as it should since the diagonal part of A commutes with the Hamiltonian. Taking a — c — 
for simplicity and requiring that the change be at least as large as the uncertainty leads to, 

I sine cose {be''*' + b*e-''t') | 1 

This simple example illustrates that a discrete evolution of the type being considered 
could be independent of /3, could have an oscillatory n-dependence and there could be a 
sub-class of states which are not eigenstates of the Hamiltonian and yet will not exhibit 
detectable evolution. These are of course special properties of the particular system. 

More generically, one could try to see the (3 dependence in the limit /9 ^ 0. The expec- 
tation values can then be expressed as a power series in (5. It turns out that, 

{A)2m ~ (A)o + o{(3^) 

(AA)L ^ i^A)l + oiP') (34) 

Thus the difference of expectation values at n = 2m -|- 1 and n — 2m — 1 is of order 
while the uncertainty at n = 2m is the uncertainty at n = plus a term of order f3'^. For 
detectability then the uncertainty at n = must be comparable to This gives a hint 
about the role of If we select a set of observables with respect to which we wish to detect 
an evolution then /3 should be chosen to be of the order of (or larger than) the square root 
of the uncertainty of the observable with the smallest uncertainty. Note that this gives a 
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lower limit on the value of (3. In this manner, the criterion of detectable evolution can be 
used to get some condition(s) on /3. However there does not seem to be a simple way of 
obtaining an analogue of uncertainty relation. For continuous family |$(t)), there is no P, 
the evolution is unitary and the usual results follow. 

Remarks: 

1) The possibility of a discrete time arises naturally for a theory presented in a frozen 
time form. Even conventionally presented theories can be cast in this form and we exploited 
this in constructing our example. For such theories (always a constrained theory), one 
needs to choose a suitable degree of freedom as a 'time degree of freedom'. One can always 
view the kinematical Hilbert space as a tensor product of Hilbert space of the time degree 
of freedom and the Hilbert space of the rest of the degrees of freedom. Solutions of the 
constraint can then be obtained as a families of vectors in the non-time sector. The 'time' 
now appears as a label for each of the family and this could be continuous or discrete. The 
families themselves are then determined as solutions of differential or difference equation. 
The form and order of the equations depends on the form of the constraint, the choice of 
representation (or choice of basis in the Hilbert space of the time degree of freedom) and of 
course on the choice of the time degree of freedom. Except for these details, a discrete time 
presentation can be set up generally. 

A continuum approximation for a discrete presentation can be looked for in the usual 
manner as indicated in subsection 11 B. The stronger notion of 'prc-classicality' however 
may not always be realizable. Even in our case, we do have asymptotic solutions which 
do have a pre-classical limit but exact solutions do not admit such a limit. Emergence of 
continuous time however can be sought via a transform instead of a limit. 

Note that we did not need to be explicit about the Hilbert space of non-time degrees 
of freedom. Even the number of degrees of freedom is unimportant for discrete time 
description. These details are of course crucial in the construction of Hq and observables. 

2) The parameter /3 is a priori completely arbitrary and has dimensions of inverse energy. 
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h(3 is thus a time scale. What possible interpretation can one ascribe to this parameter? In 
particular, does it reflect some intrinsic property of the system (and thus is selected by the 
system) or is it related to the resolutions with which a set of measurements are performed 
(and thus is to be selected by the experimenter)? Quite independently, there are two time 
scales: the intrinsic one by which a physical system keeps evolving and the clock scale eg. 
the least count of an actual clock of an observer. 

If hp is the intrinsic time scale then logic of continuum approximation would require 
that the clock scale be much much larger that the intrinsic one so that continuum time 
description is a very good approximation. Conversely if the clock scale is comparable to 
the intrinsic scale then one should use the discrete evolution. The schematic argument for 
small /? given above would now imply that discrete evolution may still not be observable if 
the uncertainties in the tracked observables are larger than permitted by the (3. 

The intrinsic time scale may be roughly estimated to be of the order of the inverse of 
the maximum uncertainty in energy measurement that may occur in the system. This 
need not be infinite, since there would be a maximum energy above which modeling of 
the physical system breaks down - eg. particle a box would not be a valid description 
for arbitrarily high energies though the spectrum of the Hamiltonian is of course unbounded. 

If however hfi is not an intrinsic time scale, then it needs to be adjusted depending 
upon what observables are used for tracking evolution. The small argument gives lower 
bounds on f3 in terms of the uncertainties of the tracked observables. Unlike the continuous 
evolution, which is also unitary, the detectability of evolution is directly dependent on 
particular observables used for tracking. 

We are unable to decide between the alternatives. It is possible that hp is a scale 
intermediate to the intrinsic and the clock scales. Considering analogy with the Immirzi 

parameter in the context of LQC, interpretation of (3 here may also throw some hints about 
the role of the Immirzi parameter. 

3) One of the powerful methods of studying semi-classical limit for systems with 
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phase space i?^^, is via the Wigner distribution function |TT|. It is essentially a dou- 
ble Fourier transform of expectation value of a certain unitary operator. In the usual 
continuous time description, the distribution function satisfies the classical Liouville 
equation to leading order in h when the expectation value is taken with respect to a 
solution of the Schrodinger equation |T^. It would be interesting to repeat the steps with 
discrete time though it looks complicated because of the structure of the difference equation. 



4) This work has been motivated by the LQC work. So what does it say about discrete 
time evolution in LQC? As mentioned in the introduction, at present the discussion of 
physical quantities such as inner products and matrix elements of observables is at a 
schematic level. The present work points out the possible pitfalls one may encounter. 
Our analysis of pre-classical limit indicates that such a limit could exist, at the level 
of states, only for solutions whose asymptotic behaviour has a characteristic root equal 
to 1. It may still go through at the level of expectation values if the largest root is 
positive. For the isotropic cosmology with positive spatial curvature (Bianchi-IX), for the 
expressions given in 0, there is neither a root equal to one nor is the highest root positive. 
However, independent of whether a pre-classical limit exist or not, one could go ahead 
with a discrete equation at the level of expectation values. One may then take recourse to 
Wigner distribution formalism to explore the semi-classical limit. This of course needs the 
Wigner distribution formalism to be developed in the context of the polymer representation. 



While one may construct families of states and even arrange schemes to distinguish cor- 
responding expectation values, it leaves unanswered the question as to why does any system 
evolve at all? 



Acknowledgments 

I would like to thank Martin Bojowald for useful comments and pointing out the earlier 
works on discrete time quantum mechanics. 



[1] M. Bojowald 2001, Class. Quant. Grav. 18, 1055; |Kr-qc/0008052 



19 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 



M. Bojowald 2001, Class. Quant. Grav. 18, 1071; |gr-qc/0008053 . 



M. Bojowald 2002, Class. Quant. Grav. 19, 2717; fer-qc/ 0202077 



M. Bojowald 2001, Phys. Rev. Lett. 86, 5227; |gr-qc/010206g . 
M. Bojowald 2001, Phys. Rev. Lett. 87, 121301; fer-qc/0104072 



M. Bojowald 2001, Class. Quant. Crav. 18, L109; |gr^/0105115 . 

M. Bojowald 2002, |gr-qc/020605^ . 

T. D. Lee 1983, Phys. Lett 122B, 217; 

R. Friedberg and T. D. Lee 1983, Nucl. Phys. B225 [FS9], 1; 

C. M. Bender, K. A. Milton, D. H. Sharp, L. M. Simmons,Jr and R. Stong 1985, Phys. Rev 
D 32, 1476; 



C. M. Bender, L. R. Mead and K. A. Milton 1993, |hep-ph/9305"24^; 



G. Jaroszkiewicz and K. Norton 1997, J. Phys. A 30, 3115, |hep-th/970307g| ; 



G. Jaroszkiewicz and K. Norton 1997, J. Phys. A 30, 3145, |hep-th/970308q ; 
G. Jaroszkiewicz and K. Norton 1998, J. Phys. A 31, 977, [hep-th/9707029| ; 



G. Jaroszkiewicz and K. Norton 1998, J. Phys. A 31, 1001, |hep-th/970703q ; 



G. Jaroszkiewicz and K. Norton 1998, [hep-th/9804165| ; 

C. Di Bartolo, R. Gamibini and J. Pullin 2002, |gr-qc/0205123 . 
[9] A. Ashtekar, R.S. Tate and C. Uggla 1993, Int. Jour. Mod. Phys. D2, 15; |gr-qc/9302027 . 
[10] Saber N. Elyadi 1996, An Introduction to Difference Equations (Springer- Verlag, New York). 
[11] M. V. Berry 1977, Phil. Trans. Roy. Soc. A 287, 237. 
[12] J. J. Halliwell 1987, Phys. Rev. D 36, 3626. 



20 



